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It is shown that a new type of the self-similar spherical ionization waves may exist in gases. All 
spatial scales and the propagation velocity of such waves increase exponentially in time. Conditions 
for existence of these waves are established, their structure is described and approximate analytical 
relationships between the principal parameters are obtained. It is notable that spherical ionization 
waves can serve as the simplest, but structurally complete and physically transparent model of 
streamer in homogeneous electric field. 



PACS numbers: 51.50.+V, 52.35.-g, 52.80.-s, 72.20.Ht 

The impulse electrical breakdown of various media is 
usually caused by appearance and propagation of the ion- 
ization waves (IW), in particular, the streamers. How- 
ever, the theory of this important and interesting subject 
of long-term investigations is still far from the comple- 
tion. The reason for this state of affairs is quite obvi- 
ous: it is necessary to solve multidimensional and non- 
linear system of balance equations together with Pois- 
son's equation for description of IW even in the frame- 
work of simplest hydrodynamic approximation. So the 
analytical description is possible only if they are highly 
symmetrical objects that possess the property of self- 
similarity. Among them stationary plane IW are the sim- 
plest ones and some rigorous results have been obtained 
for them [l|, 0, S] ■ Meanwhile the most of IW, studied 
experimentally and by numerical simulation, were nei- 
ther stationary nor plane. Therefore subsequent search 
for the self-similar solutions, describing non-stationary 
curved IW, is of extreme importance. In this field first 
results were obtained recently: the authors of papers Q 
studied cylindrical and spherical IW in gases, assuming 
the absence of full current and, therefore, the constancy 
of the wave charge value Q. This condition causes the 
reduction of the maximum field strength Em = \E(rt)\ 
at the front with rise of its radius rf in conformity with 
the same power law as for the initial undisturbed field. 
Such IW slow down and attenuate in the course of time 
and their evolution becomes self-similar only at so small 
Em, that the impact ionization becomes negligible and 
the wave propagates only due to the electron drift. 

In the present work I will examine a modified problem 
by assuming, in contrast to A], that the current flows 
through the wave and supports the constancy of the sur- 
face charge density Q/Anr 2 of the spherical wave front 
and, hence, Em- It will be shown, that in this case 
there exist asymptotically self-similar solutions, which 
describe undamped and exponentially accelerating IW. 
Such waves can serve the simplest, but structurally com- 
plete and physically transparent model of streamers in 
homogeneous electric field, for which the constancy of 
Em is the most plausible hypothesis 

Let us consider IW which propagates inside the spheri- 



cal gas-filled capacitor towards the external electrode (for 
the certainty - to the anode). In the framework of "mini- 
mal" hydrodynamic model @, Q the distributions of elec- 
trons nit, r) and positively charged ions p (t, r) within the 
wave are described by the system of balance equations 
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where the operator V r = r gf' ■ The impact ioniza- 
tion coefficient a, drift velocity v and diffusion coefficient 
D of electrons are considered as assigned instantaneous 
and local functions of the field strength E. Usage of the 
first of equalities (p} and Poisson's equation give the con- 
servation law of the full current J (t) along the axis r 
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where q is the elementary charge, e is the dielectric con- 
stant of the medium. Our goal is the search for solutions 
of the form /j (x) fi (t) > where self-similar coordinate 
X = ip (t/to) r/ro, to and To are the positive constants 
and index i — n,p, E. It is possible to show, that allow- 
ing for impact ionization with any plausible dependence 
a (E) this task is solvable only if if; (t/to) = exp (—t/to) 
and every ifi (t) = 1. In this case all addends in equations 
(HJ), (J2J), that contain operator V r oc exp (—t/to) V x , de- 
crease exponentially in the course of time and at t 3> to 
our system takes the asymptotically self-similar form 
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Here the following dimensionless variables and parame- 
ters are introduced: time 9 = avt, concentration a (x) — 
qn/eaE, field strength F (x) = —E/E, drift velocity 
V (F) — v (E) /v, impact ionization coefficient A(F) = 
a (E) /a, t — avto, v — v(E), ro — r c (jo/ 'eavE) 1 / 2 , 
where a and E are constants, which characterize the de- 
pendence a (E) (see further), r c is the cathode radius, jo 
is the initial current density at the cathode. During the 
derivation of equations ([3|) , (J4j) it has been assumed that 
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J(t) — —4:TTrlj exp(2t/to), since condition fE{t) = 1 
is satisfied ahead of the front only if waves charge (and 
hence cathode current) grows according to the law Q oc 
exp(2£/io). The number of equations was reduced by one 
in comparison with input system, since the space charge 
density q(p — n) oc exp(— t/to)V x F\ t ^oo — * and there- 
fore n rs p for all X- Physical reason for this lies in the fact 
that at large t the front speed Uf — rf/to oc exp(t/to) 
becomes much greater than average electrons drift veloc- 
ity. For the same reason equations ((3]), (j4|) turn to be 
invariant with regard to simultaneous sign reversal of E 
and jo- Consequently the parameters of self-similar IW 
do not depend on the electric field direction ||. 

The right sides of equations ([3]), (j4|) and their second 
derivatives are continuous in terms of the variables x> u i F 
at x > and F > 0. Consequently solution of these equa- 
tions exists and it is uniquely determined by boundary 
conditions [6| . I will suppose that the radius of external 
electrode (anode) r a 3> rj and that there exist back- 
ground charge carriers with the concentrations n + = p + 
ahead of the front. In this case the boundary conditions 
can be written as 



<t(oo) = er + , = 0. 



(5) 



The complete solution of the boundary-value problem 
©-([5]) can be obtained only by numerical methods. To 
attain this goal I applied the Runge-Kut method of the 
4th order with the adaptive grid, using the indepen- 
dent variable y — 1/x and solving the problem in the 
finite segment < y < 1/Xc- The right boundary of 
this segment 1/Xc = (jo/z&vE) 1 / 2 exp(6*/r) corresponds 
to value r = r c , i.e., the surface of internal electrode. 
Calculations were carried out using the simplest empir- 
ical approximations v = fx \E\ (p, is the electrons mobil- 
ity) and a — a exp(— E/\E\), for which V(F) — F and 
A(F) = cxp(— 1/F). The use of any other physically 
justified dependencies v(E) and a(E) gives only quanti- 
tative changes of results. 

Typical distributions of field and concentration in 
spherical IW are presented in Fig. |TJ Results for other in- 
stants of time can be obtained by a simple shift along the 
logarithmic x-axis. It is evident, that the wave consists 
of four qualitatively different regions f2 c , f2_, £lo and f2+. 
In the region f2 + ahead of the front the conductivity is 
negligible, therefore it follows from ^ that 
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This formula describes field distribution around the 
charged sphere. Substituting ([6]) in ([3]) one can get after 
integration the interrelation between concentration and 
field strength 



a w cr+ exp | ~ [770 - 77 (F) } j , 



(7) 



where 77 (F) — Jp M AdF and ryo — V (0). These formulae 




FIG. 1: Structure of spherical ionization wave at the moment 
9 = 5r with jo = 5 • W~ 5 eaEH, r = 300 and a+ = 10" 6 . 



are valid until err <C 2. Fast growth of conductivity at 
the wave front leads to the appearance of the region Qq, 
in which the field first reaches maximum value Fm with 
X = Xm and then begin decrease. In this region both 
suppression of field due to the plasma conductivity and 
abrupt change of this conductivity due to intense impact 
ionization are important factors. The thickness S of re- 
gion Qq usually much less then xm and the field strength 
changes insignificantly in flo. Therefore it is possible to 
assume x = Xm and V(F) = Fm in equations (J3|), ((!]), 
but sharp exponential dependence A(F) should be al- 
lowed. For this case the solution of (J3j , (J4j) was obtained 
in paper [3| and with V(F) — F it take the form 
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where <jm — & (xm)> and A is positive (with x > Xm) or 
negative (with x < Xm) root of the equation 



A - 1 + exp(-A) = r}j(JM- 



(10) 



Field reduction behind the wave front causes the forma- 
tion of the region fl _, where impact ionization is negli- 
gible and the solution takes the form 
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where F m — F (Xm) is minimal field in the wave. With 
X < Xm the field strength increases in the direction to 
the cathode, that is caused by the need for conservation 
of the full current through plasma with constant conduc- 
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tivity a = <7_ and with decreasing area 4irr 2 oc x 2 ■ Since 
usually to — 3> 2, the field increases as F w F m (xm/x) 2 
until in the region f2 c (% < recurring growth of a 
begins due to the impact ionization. As a result of dras- 
tic dependence A (F) sufficient growth of a is ensured 
due to the small increase of F. Therefore near the cath- 
ode field strength approaches to constant value F c , while 
conductivity increases according to power law: 

F«F C , a « a c ( X ) = (13) 

If F c is the root of equation tF c A(F c ) = 2, then formulae 
(fTB")) give one of exact solutions of the equations (O, Q. 
This solution does not satisfy to boundary conditions {5} , 
but it is attractor for the real solutions with \ ~~ * 
(see Fig. [TJ. On the cathode itself the conductivity 
a = l/(F c Xc) and increases exponentially with the incre- 
ment 2 /to in the course of time, as it must be to provide 
the same increase of the current. 

In order to obtain complete description of IW it is nec- 
essary to join the solutions (f6"))- (fl"3")) between themselves. 
Taking into account equality (JmFuXm = 1 the formula 
© can be rewritten in the form t<jm — 2(\om I Xm) 2 , 
where Xom — y t /2Fm by definition. Smallness S makes 
it possible to disregard a difference between (xom / Xm) 2 
and 1. As a result the estimations of concentration 
(Xm ~ 2/t, front position Xf = Xm ~ \Jt/2F m and 
front speed Uf — rf/to = (roto^ 1 y / T~/2FM^ exp(t/t ) 

are derived immediately. Joining dependencies j7|), ((9|) 
and taking into account the fact that ahead of the front 
A ~ (ry/(TM + l) ~ (T-77/2 + 1), it is easy to obtain approx- 
imate equation 

1 ' ■ (14) 
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that establish the relation of maximum field Fm with 
control the parameters er+, r. Formula for the electron 
concentration behind the front <r_ follows from © : 



■ exp (—A- 



(15) 



where A_ is the negative root of equation 
[A_ - 1 + exp (-A_)] = rjoT/2. Instead of (15]) it 
is possible to use the simple approximation 

9 

a- = a a rjo, a a = 1 + — (16) 

2 In 770/0 + 

that differs from (fT"5)) less than to 2% under condition 
rjo > 15cr-|-, which is satisfied well usually. The analysis 
of the numerical solutions of the boundary- value problem 
©-© shows, that (xqm/xm) 2 < 1.1 with a + < 10" 4 . 
This inequality justifies the approximations xom = Xm 
used above. The stright comparison with numerical re- 
sults shows that formulae (fT4 |) -(fT6 |) are carried out with 
a good accuracy in all examined cases, as it is evident 
from Fig. 2,3. Further, equalizing conductivity currents 
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FIG. 2: Plot of maximum field strength against response time 
with different values of a+. Points: numerical solution, line: 
calculation according to (Tl4)) . 




FIG. 3: Plot of electron concentration behind the front 
against maximum field strength with different values of a+. 
Points: numerical solution, lines: calculation according to 

(Hsu. 



on the surfaces x — Xm and x — Xm, one can obtain 
relationship for estimation of F m : 



<r-F m = 2a x F M /r, 



(17) 



where a\ 



(Xom I Xm) 2 - Formula (|17p ensures high accu- 
racy if a x = 1.15acr (see Fig. 4). Finally, join of solutions 
(fTT|l - (fT3|) give the position of boundary X- ~ (o'-F'c) -1 / 2 
between regions tt c and f2_. 

The obtained formulae establish interrelations of all 
principal values, that define the structure of IW (F c , x~, 
Xm, cr-, F m , Xf, Fm), with the control parameters (0 + , 
t). Thus they give the complete description of waves 
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FIG. 4: Plot of conductivity current density behind the front 
against displacement current density ahead of the front with 
different values of cr+. Points: numerical solution, line: cal- 
culation according to (|17[) with a x — 1.15a CT 

for large t. It is clear, however, that the acceleration of 
IW cannot continue infinitely for a long time. Firstly, 
waves must be impeded due to the appearance of vor- 
tex electrical field, when the front velocity approaches to 
the velocity of light [3]. Secondly, exponential build-up 
of plasma conductivity near the cathode sooner or later 
will be suppressed by nonlinear recombination, electron- 
ion scattering, saturation of cathode emission ability or 
due to the total ionization of medium. The above men- 
tioned and similar mechanisms will fail to appear, if a 
radius r a of external electrode is not too great. In this 
case self-similar propagation of the wave will continue 
up to the front contact with the external electrode at 
the time t « to \n(r a /r Xf), if external electric circuit 
ensures the exponential build-up of current. It is very 
difficult to attain this in practice since value r = 1 corre- 
sponds to response time to of the order of several picosec- 
onds. Therefore a reasonable question arises concerning 
the existence of the real objects, for which a model of 
self-similar accelerating IW is descriptive. 

In our opinion such objects are, in particular, the 
streamers in uniform electric field, which repeatedly were 
studied by two-dimensional numerical simulation. The 
most detailed data are represented in the papers Q for 
nitrogen at normal pressure. Analysis of these data indi- 
cates the qualitative coincidence of structure and evolu- 
tion characteristics of streamer and spherical self-similar 
IW. This coincidence is so complete that (being dis- 
tracted from the nicety of the streamers transverse struc- 
ture) it is possible to declare: streamer in uniform elec- 
tric field is simply a cone with spatial angle at vertex of 
order 0.17T, that have been "cut out" from self-similar ac- 
celerating spherical IW. However, the evolution of such 
streamers inside a large parallel-plate capacitor possesses 



a number of features in comparison with spherical IW. 
First of all, an exponential growth of the charge with 
the increment 2 /to (which is necessary for the self-similar 
propagation) is ensured not due to external control, but 
automatically during evolution of streamer as the self- 
organizing object, that obtains energy from the uniform 
field. Secondly, in gases the seeding electrons ahead of 
the streamers front arises mainly due to ionizing radi- 
ation, emitted from plasma behind the front. In this 
case IW can propagate almost self-similarly as long as 
rfctph < 1, where a p h is minimal absorption coefficient 
of the ionizing radiation. This inequality was fulfilled in 
Q, where a value a p h = 5 cm -1 was used in calculations 
up to the value r/ ~ 0.2 cm. Thirdly, streamers evolu- 
tion must lead to an increase in the current in external 
circuit exponentially with the increment 3/io, although 
the increment of charge growth is equal to 2 /to just as in 
the spherical wave. This nontrivial result was revealed, 
but not explained in [7fl . It follows directly from Sato the- 
orem Q and from the condition for exponential growth 
of all streamers spatial scales with the increment 1/to- 

In conclusion it should be noted that self-similar IW 
can also exist even if uniform distributed seeding charge 
carriers are absent (i.e., at a + = 0), but some addi- 
tional mechanism of medium ionization exists ahead of 
the front. The effective action radius r g of this mecha- 
nism, obviously, must increase in proportion to r^. This 
scenario can be realized in semiconductors, where tunnel 
ionization is the basic mechanism of electron-hole pairs 
generation ahead of the front Q . The rate of this process 
depends on E locally. Hence condition r g « rj is satis- 
fied automatically and streamer in the uniform field must 
propagate in accordance with the described scenario for 
a reasonably long time. 
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